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AN ANALOG STUDY OF THE RELATIVE IMPORTANCE OF 
VARIOUS FACTORS AFFECTING ROLL COUPLING 
By Joseph Weil and Richard E. Day 


SUMMARY 


An analog study of the roll-coupling problem has been made for a 
representative swept -wing and tailless delta-wing configuration. The 
investigation, made primarily for subsonic flight conditions, included 
the determination of the effects of wide variations in many of the per- 
tinent aerodynamic derivatives on the motions developed in rolling 
maneuvers. The influence of large changes in principal axis inclination 
an d mass distribution is also considered. 

The results indicated that as first predicted in NACA TN 1627 the 
relationship between the longitudinal and directional stability is of 
paramount importance. For most current designs an optimum condition 
exists when the natural frequencies in pitch and yaw are approximately 
equal. Increases in pitch damping had a pronounced favorable effect in 
reducing the amplitudes of the motions encountered and were, in general, 
considerably more effective than corresponding increases in yaw damping. 

Practical redistribution of mass produced only relatively minor 
changes in the overall results . 

It was found that the amplitude of the motions developed for a given 
aileron deflection depends to a large extent on the duration of the 
maneuver (change in bank angle). Limited studies indicated that 90 ° roll 
maneuvers would be considerably less critical than 360° rolls. The angle 
of attack of the principal axis has an important bearing on the behavior, 
particularly in the absence of other disturbing functions. If the ini- 
tial angle of attack is maintained constant, a reduction in altitude will 
delay critical conditions to a higher roll rate but the maximum ampli- 
tudes may be only slightly affected. 

Small inadvertent stabilizer inputs can greatly affect the magni- 
tude of the motions developed. 

It is difficult to generalize on the effects of Mach number varia- 
tion because this variable affects many of the controlling parameters. 
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Utilizing simple concepts proved useful in assessing the qualitative 
effects of many of the aerodynamic and inertia parameters and changes in 
flight condition. The calculated lover resonant frequency generally cor- 
responded to the average roll velocity at which the more serious motions 
could he expected. 

It is indicated that rational design procedure can avoid the problem 
of serious roll coupling at supersonic speeds and minimize the problem 
in the subsonic speed range for the configurations of the type considered. 


INTRODUCTION 


Until recently the motions of a rigid airplane with deflected ailerons 
could be adequately predicted in most instances from the classic three- 
degree -of -lateral freedom solutions (for example, ref. l). That the 
longitudinal modes can exert a powerful effect on the overall motions of 
the rolling aircraft was first demonstrated theoretically in 1948 (refs. 2 
and 3)j however, not until 1954 were strong manifestations of coupled 
motions in rolling maneuvers experienced on a full-scale airplane (ref. 4). 

Following the flight experiences described in reference 4, a five- 
degree-of -freedom analog study was initiated by staff members of the 
NACA Langley Laboratory and NACA High-Speed Flight Station to determine 
whether the actual flight motions could be calculated and also to deter- 
mine the effects of variations in certain of the stability parameters. 

Some of the results of these prelim i nary studies were reported briefly 
in references 5 and 6. 

An analog investigation of several generalized airplane configura- 
tions encompassing a considerably broader scope than the work of refer- 
ences 5 and. 6 has recently been completed at the NACA High-Speed Flight 
Station. In this study wide variations in many of the pertinent aerody- 
namic derivatives were investigated at subsonic and supersonic speeds. 

The effects of large changes in principal axis inclination and mass dis- 
tribution are also included. 

The primary purpose of this paper is to summarize the information 
obtained from the analog calculations and to compare the results with 
the trends predicted from a slightly modified version of reference 2. 

SYMBOLS 
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wing span, ft 
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c m 


Cn 
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h P 

it 

% 

I X e 

x xz 

% 

x z 

M 

M 

m 

N 

P 

P 


lift coefficient, y 

^PV 2 S 

rolling-moment coefficient, moment 

|pV 2 Sb 

pitching -moment coefficient, Pitching moment 

■ipV 2 Sc 

yawing-moment coefficient, having moment 

|pV 2 Sb 

lateral-force coefficient, La_t e ^ a l force 

|pv 2 s 

wing mean aerodynamic chord, ft 
acceleration due to gravity, ft/sec 2 
pressure altitude, ft 
stabilizer deflection, deg 

moment of inertia of airplane about X-axis, slug-ft 2 

moment of inertia of rotating engine parts about X-axis, 
slug -ft 2 

product of inertia referred to X- and Z-axes, slug-ft 2 
moment of inertia of airplane about Y-axis, slug-ft2 

moment of inertia of airplane about Z-axis, slug-ft 2 
Mach number 

pitching moment, ft-lb 

airplane mass, — , slugs 
8 

yawing moment, ft-lb 
rolling velocity, radians/sec 
average rolling velocity, radians/sec 
steady rolling velocity, radians/sec 
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q. 

r 

S 

t 

V 

w 

X, Y, Z 
a 

% 

a 0 


ap 

p 

Ax, AP 
5 H 




P 


pitching velocity, radians/sec 
yawing velocity, radians/sec 
wing area, sq ft 
time, sec 

true airspeed, ft/sec 
airplane weight, lb 
body axes of airplane 

angle of attack of airplane body axis, radians or deg 
angle of attack at zero lift, radians or deg 

initial angle of attack of airplane body axis, deg 

initial angle of attack of principal axis, deg 
angle of sideslip, radians or deg 

increments from initial conditions, radians or deg 

total aileron deflection (positive for right rolls), radians 
or deg 

rudder deflection, radians or deg 

angle between body axis and principal X-axis, positive when 
reference axis is above principal axis at the nose, deg 

fraction of critical damping in pitch of nonrolling aircraft, 

-cMq 

AW-McJy 

fraction of critical damping in yaw of nonrolling aircraft, 
-bN r 

Wn^T 

mass density of air, slugs/cu ft 
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Ap 

CD 0 


increment in angle of bank, deg 

rotational velocity of engine rotor, radians/sec 

nondimensional undamped natural frequency in pitch of non- 
rolling aircraft (ratio of pitching frequency to steady 
rolling velocity) 

nondimensional undamped natural frequency in yaw of nonrolling 
aircraft 


Ct , Ct , , Ci , Ci. , 

■MO. L it z 5 at 

Cmo,? C m p, C m g e , Cm^, 

> indicates derivative with respect to subscript 

c np> c ng > Cn 6r' 
c Yp' Ma' 

C V c Zr , c np , c nr , 

V C Y r > N r 


indicates derivative with respect to 
2^ X subscript 


c mq> Mq 


indicates derivative with respect to 
X subscript 


Dot over a symbol indicates derivative with respect to time. 


CALCULATIONS 
Motion Studies 


The basic time histories, upon which the analysis of this report is 
based, were calculated in real time using a five-degree-of-freedom approach 
(forward speed assumed constant). The fundamental equations used to cal- 
culate the motions are shown in table I. These equations were solved by 
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means of a Goodyear Electronic Differential Analyzer (GEDA). 1 The initial 
conditions used in the calculations simulated either trimmed level flight 
or constant g accelerated flight. Aileron deflection was generally the 
sole control input used in the calculations. 


Several thousand time histories were run off in the course of the 
general study. The initial calculations were made with the simple aileron 
input of figure 1(a) (input A). A ramp-type input of 50° per second was 
used until the desired control angle was reached. At the nominal hank 
angle (usually 360°) the control was returned to zero at 50° per second. 

It was found in a number of instances that this simple aileron input pro- 
duced unrealistic results (for example, fig. 1(b)) because it failed to 
arrest the rolling motion after the aileron was returned to neutral. To 
obtain a more realistic evaluation in these instances it was necessary 
to use a control stick in conjunction with the analog, as shown in 
figure 2. A typical control input using the stick (input B) is shown 
in figure l(c). It should be noted that the inputs A and B were similar 
until the point of corrective control. In the latter instance, the oper- 
ator attempted to stop the rolling motion as rapidly as possible. 

A number of approaches were tried in an effort to arrive at the most 
s iSfl-ificant parameters for summary purposes. It was finally decided that 
the maximum positive and negative excursions in angle of attack and side- 
slip (fig. 1(c)) plotted against the average roll rate obtained in a given 
maneuver might provide the best overall representation. The average roll 
velocity is defined as the bank angle at control reversal divided by the 
time required to reach this angle. The bank angle was determined from 
the direction cosine m^ (table I). 


Simplified Analysis 


Throughout the analog program reference 2 served as a valuable guide 
notwithstanding the fact that the a na lysis of reference 2 assumed constant 
roll rate and made other simplifying assumptions. It was thought desir- 
able, however, to account for engine gyroscopic effects and to use the 

actual value for the ratio — — X instead of assuming uni t.y. 

X Y 


Another approach in considering the effects of engine momentum on 
rolling maneuvers is presented in reference 7. 

It is shown in appendix A that inclusion of the previously mentioned 
modifications to the theory of reference 2 results in an increment to the 

^The GEDA equipment was made available through the cooperation of 
the Air Force Flight Test Center. 
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parameters used to define the stability charts and also lowers the hori- 
zontal boundary slightly (appendix A and fig. 3)* The basic meaning of 
figure 3 Is identical to that discussed in reference 2 in that the region 
between a boundary and axis indicates a divergent condition. The roll 
velocity at which the effective stability first becomes unstable for a 
condition of zero damping is referred to as the resonant frequency and 
is given by the following equations 


C np | pV 2 Sb I Xe <% 
^ ± ]-TT^r + 2~(£ Y ~'lx] 


( 1 ) 


applicable where the vertical boundary is first intersected as the roll 
velocity is progressively increased, and 


■*1 


-C 


% 


\ pV 2 Sc 




L x 


2 ( I z - *x) 


( 2 ) 


applicable where the horizontal boundary is first intersected. 

Note that equations (l) and ( 2 ) are identical to equations 10 and 
11 of appendix A} however, the notation p Q of reference 2 was retained 
in the appendix whereas the average roll velocity is used in the remainder 
of this paper. 


SCOPE 


Two basic fighter airplane configurations are included in this study. 
One configuration had the mass and aerodynamic characteristics represen- 
tative of current swept-wing airplane types, while the other used param- 
eters estimated for a tailless delta-wing type. All the rolling maneuvers 
for the swept-wing airplane included in this paper were calculated for 
a Mach number of 0.7, whereas the delta configuration is presented for 
M = 0.8 and M = 1.2. Unless otherwise indicated, the basic aerodynamic 
and mass parameters shown in table II were used in the calculations. The 
aerodynamic derivatives were based on wind-tunnel results, flight tests, 
and theory. The derivatives used (with the exception of Cj ) were assumed 


invariant with angle of attack. Table II presents the variation of C’ 

L 

with a used for the several basic flight conditions. 
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For each configuration a series of calculations were made with suc- 
cessive increments in aileron deflection using the basic derivatives of 
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table II. Additional series were obtained by varying the stability deriva- 
tives and C n(3 and the damping derivatives C m ^ and over 

a rather wide range. The effects of roll direction and duration and 
initial angle of attack were also considered for each basic condition. 


To obtain a better insight into the mechanism of roll coupling, 
additional information was obtained for the swept -wing configuration with 
arbitrary variations in yawing moment due to rolling Cn.p-> dihedral effect 


C 7 , principal axis inclination, altitude, and mass distribution. The 

importance of relatively small inadvertent stabilizer inputs was also 
determined. 


PRESENTATION OF RESULTS 


A summary of the figures presenting the results of this investigation 

are: 


Effects of: 

Cm at various levels of C n „ 

“bo p 

Roll direction 

Roll duration 

Damping in yaw * 

Damping in pitch 

Inadvertent stabilizer input 

Initial angle of attack and 

principal axis inclination 

Altitude 

Yaw due to roll ' 

Dihedral effect 

Mass distribution 

Analysis of coupling problem on 

specific design 


Figures 

4 to l6 

17 and 18 
19 to 21 
22 to 25 
26 to 33 
34 and 35 

36 to 40 
4 l 

42 and 43 
44 

45 to 47 
48 


DISCUSSION 
Individual Effects 


Importance of C^, C n(J relationship.- Figures 4 to 6 present a 
summary of the effects of variations in C^ at three levels of C n p 
for the swept-wing configuration. The data of figures 4 to 6 present 
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the envelopes of the maximum angle of sideslip and angle of attack excur- 
sions from trim (fig. 1 (c)) as a function of average roll velocity in 
nominal 360° left rolling maneuvers. The results indicate that the value 
of Cidq, has a profound effect on the calculated amplitudes at all levels 
of directional stability. Moreover, the value of at which the 

smallest angles are obtained appears to increase as C n p is increased. 


The average roll velocity at which the peak amplitudes occur would appear 
to vary in a somewhat erratic fashion. To analyze the trends shown, the 
simple theory of reference 2 (modified as indicated in appendix A) was 
utilized. In figure 7 the nondimens ional frequency parameters in pitch 
and yaw are plotted with suitable terms added to account for engine gyro- 
scopic effects. For any given flight condition, variation in the roll 
velocity will trace out a curve which for the values of the aerodynamics 
and engine momentum considered is very close to a straight line. For 
conditions falling between the coordinate axes and the vertical neutral 
stability boundary (region A) the simple theory predicts primarily a 
directional divergence, whereas conditions between the coordinate axis 
and the horizontal neutral stability boundary (region B) are primarily 
indicative of a pitch divergence. It should be noted that, inasmuch as 
the physical and inertia characteristics were assumed constant, the slopes 
of the lines shown in figure 7 are for all practical purposes dependent 
on the ratio of to C n p. Each condition in figures k to 6 is num- 

bered and represented in figure 7bya line. The several instances of 

identical C m /C ratios are indicated by a common reference number. 

%/ u p 

It can now be seen that the farther from the origin the representative 
line intersects the boundaries of figure 7> the worse the coupled motions. 
Thus, conditions 1 and 5 which represent the most extreme ratios of natural 
frequency in yaw and pitch also are characterized by the least desirable 
motions. The average roll velocity at which the lower undamped resonant 
frequency for the nonrolling airplane occurs (equation 1 or 2) is shown 
ticked on the curves in figures 4 to 6. It seems that this frequency in 
general occurs fairly close to the roll velocity at which maximum excur- 
sions occur. For a given ratio of to C n ^ the higher the roll 


velocity for resonance, the greater the maximum amplitude of the motions 
for the more extreme conditions 2 and 4. Condition 3> which has a more 
desirable proportioning of stability, does not show this trend to any 
degree . 


is also evident in 
figures 4 to 6 and the results follow the foregoing analysis. 


The effect of increasing C n ^ at constant C m ^ 


Representative time histories for aileron deflections at which near 
peak amplitudes were obtained are shown for each combination of 
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and C n p in figures 8 to 10. In nearly all instances the maximum excur- 
sions occur during the recovery phase of the maneuver. Although the 
recovery is perhaps not mathematically definable, the results of fig- 
ures 4 to 6 were not particularly sensitive to the exact control manip- 
ulations used during recovery, and repeat runs produced very little 
scatter on the 'summary plots. For conditions 1 and 2 which produced an 
intersection in the upper portion of the diagram of figure 7 the predom- 
inant angle-of-attack change during the constant aileron portion of the 
maneuver was negative, whereas for conditions 4 and 5 the predominant 
angle-of-attack change was positive. The initial sideslip motion was 
always negative, although the yawing velocity histories varied markedly. 


A series of runs was made in which individual terms were eliminated 
singly from the equations of table I. The results are shown in figures 11 
and 12 for the two most extreme conditions, 1 and 5 (figs. 4 and 6 ). 


For condition 1, which represents a condition having a C n p value 
which is much too low in relation to C^, it is evident that elimination 

of the pq term in the r equation effectively de-coUples the motion 
(fig. 11(c)). When the pr term in the 4 equation is removed, however, 
little effect is shown. The elimination of the p (3 term in the & 
equation shows an effect similar to omission of the pq term. 

For condition 5> which represents a condition having a C Q p value 
much too high relative to Cn^, removal of the pr term has the primary 

de-coupling effect (fig. 12). It is seen, moreover, that the removal of 
the pq or p 3 terms produces an unfavorable effect for condition 5 * 

The importance of the C^, C n p relationship for the delta-wing 

configuration at M = 0.8 and M = 1.2 is illustrated in figures 13 
and l4. Stability diagrams for these conditions are presented as fig- 
ures 15 and l 6 . From a close inspection of figures 13 to l 6 it is seen 
that the major points brought out in the discussion of the swept -wing 
configuration are substantiated for the delta-wing configuration. 

For the M = 0.8 condition (figs. 13 and 15) the expected favorable 
effect of reducing Cn^ from -O .36 to -0.l8 (condition 1 to condition 3) 

is evident, although the magnitude of the motions is still quite large 
for the latter case. When the value of C]^ is further reduced to -0.0$), 

a considerable reduction in the peak amplitudes is obtained probably 
because of the decrease in lower resonant frequency. 


At M = 1.2, doubling the basic value of Cnp resulted in an appre- 
ciable improvement because of a more desirable proportioning of stability 
(conditions 1 and 3 > fig* l 6 ) and because of the increased stability level. 
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Detailed comparisons of the several configurations discussed should 
be made only in the light of the different initial conditions present 
(table II). The effects of some of these variables will be treated in 
later sections of this paper. 

Effect of roll direction .- The effect of roll direction is su m marized 
in figure 17 for the swept-wing configuration at C n p levels of 0.057 

and 0.1l4 per radian. The data indicate that in left rolls with the 
lower C n p> the peat motions occur at a somewhat lower average roll rate 

and are 50 percent larger than corresponding right rolls. For 

of 0.11^- a similar trend is indicated but is less marked. 

The simple theory of reference 2 modified to include the engine 
terms has been applied to the conditions of figure 17 • The stability 
chart (fig. 18) shows that for the lower C n p the representative line 

for the left roll intersects in a slightly less desirable location than 
the corresponding line for the right roll. For the higher C n p the two 

lines are coincident. The lower resonant frequencies for the undamped 
conditions are indicated by the ticks (fig. 17) and it is apparent that 
the increase in resonant frequency calculated for the right rolls cor- 
relates fairly well with the motion studies . 

It should be mentioned that although the left rolls of figure 17 
are seen to be somewhat more critical than the corresponding right rolls , 
a number of factors could cause the reverse to be true. Among these 
factors would be direction of engine rotation, and the values of Cm- , 

°at 

C n , and perhaps initial angle of attack. For all conditions considered 

in this paper however, the left roll is the more critical and all of the 
remaining results are presented for left rolls . 

Effect of roll duration .- A summary of the effects of roll duration 
is shown for the swept-wing configuration in figure 19. It is seen that 
the maximum amplitudes calculated for the 90° and l80° maneuvers are much 
lower than for the 360° rolls. This is not surprising if the motions are 
considered to be the manifestation of an effective reduction in stability 
or actual instability, thus the time duration (change in bank angle) would 
be expected to be a determining factor in the motion build up. Represent- 
ative time histories are shown in figure 20 for 90°, l8o°, 360°, and lk4o° 
maneuvers. The possible effects of roll duration sire clearly evident 
particularly for the l^ko 0 maneuver in which angles of attack and angles 
of sideslip of very large magnitude are obtained. In most instances, 
however, the amplitudes attained in the 36O 0 rolls were comparable with 
those calculated for l4ko° rolls, particularly when the recovery phase 
is included. 

According to simple theory (appendix B) it can be shown that the 
rate of divergence for an unstable condition will be directly proportioned 
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to the roll velocity and also he adversely affected when the intersection 
of the determining line with the stability boundary occurs at large dis- 
tances from the M neck n of the diagram. Thus the magnitude of the motions 
might be expected to build up quite rapidly for the delta-wing configu- 
ration (condition 1, fig. l6) at supersonic speeds. Time histories 
showing the effect of bank angle are presented in figure 21 for the delta- 
wing configuration at M = 1.2. A 3^0° roll is seen to produce a violent 
maneuver with sideslip angles greater than 20° and large a excursions. 
The l80° maneuver still develops 17-5° of sideslip and ±6° in Ax (±2g). 
Even the 90° maneuver, for this condition, is accompanied by 12° of side- 
slip but the a amplitudes are considerably reduced. 

Effect of damping in yaw .- The influence of changes in C n is 
summarized in figure 22 for the swept-wing configuration for a Cn^ 
level of 0 . 057 . A tenfold increase in C nr reduces the peak p motions 

appreciably but has very little effect on the maximum angle-of-attack 
amplitudes. A value of C nr 100 times the basic value produced no sig- 
nificant additional decrease in the maximum p amplitudes although a 
reduction in the a motion is evident. Moreover, particularly for the 
highest damping ratio, the peak amplitudes were materially increased at 
the lower roll rates. A stability diagram of the type used previously 
is shown as figure 23* From the simple theory it might be ejected that 
a more significant improvement would exist for the two conditions of 
increased damping in yaw. Note also that the peak amplitudes of figure 22 
occur near the same average roll rate, whereas the results of figure 23 
might lead one to believe the maximum excursions with increased damping 
would occur at appreciably higher roll rates. 

Recalling the results of figure ll(b), the elimination of the rp 
term in the q equation did not alter the basic level of the motion for 
conditions intersecting the vertical divergence boundary. For conditions 
intersecting the horizontal divergence boundary (fig. 12(b)) elimination 
of the rp term was very effective. 

Figures 24 and 25 show the effects of increasing by factors 

of 10 and 100 for the extreme conditions 1 and 5 (figs. 4 and 6). The 
results indicate that the increases in Cn r have a more significant 

effect in reducing the motions for condition 5 than for condition 1. 

Effect of damping in pitch .- The effects of increasing the basic 
damping in pitch for the swept-wing configuration are summarized in fig- 
ure 26 for two levels of Cnp* It is seen that a threefold increase in 
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pitch damping^ = 0 . 27 ^ results in an appreciable reduction in the 

0 motions and some reduction in the a amplitudes. A tenfold increase 
in Cmq = O. 92 ) produces a marked improvement at each stability- 

level. The stability charts of figure 27 show the effect on the diver- 
gence boundaries of increasing the pitch damping. These results correlate 
fairly well with the trends presented in figure 26 . The results of the 
time histories of figures 28 and 29 show the favorable effects of increases 
in C m q of 10 and 100 times the basic value on the motions for the two 

extreme basic conditions examined in the preceding paragraph. 

In figure 30 the effects of variations in pitch damping are pre- 
sented for the delta-wing configuration at a Mach number of 0.8. At 
this Mach number increasing the level of showed a considerable 

improvement but was not quite as effective as for the swept -wing con- 
figuration. The stability diagram is presented in figure 31 for this 
flight condition. 

At a Mach number of 1.2 the results shown in figure 32 indicate a 
relatively small effect of increasing Cm^ until an extremely large 

level is assumed. It is believed this condition is caused by the inter- 
section of the basic determinant line with the divergence boundary at a 
relatively large value of the ordinate. It is seen (fig. 33) that for 
this condition moderate increases in Cm^ would be expected to have 

only a slight favorable effect. 

Effect of stabilizer input .- In the previous section the important 
effects of pitch damping were discussed. Inasmuch as these changes are 
attributable to reductions in the pitching velocity, it might be reasoned 
that even small stabilizer inputs during the rolls could likewise have 
an important effect on the results. Figure presents a summary of 
results for simultaneous stabilizer and aileron inputs for C n p = 0. Il4. 

The type of stabilizer input used is illustrated in the representative 
time histories of figure 35* It is evident that 1° of airplane nose-up 
stabilizer results in a 50 -percent increase in the amplitudes of the 
0 motions and also produces large increases in the a excursions, 
whereas a stabilizer input of the opposite sense has a large alleviating 


%he damping ratio considering only C m q is shown in the figures 

because this value was used in calculating the boundaries shown. Of 
course, the total angle of attack damping ratio is also dependent on Cn^ 

and the lift-curve slope. Thus, for the swept wing airplane (figs. 26 

and 27 ) the total damping ratio will be about 0.l6 higher than £q, whereas 

for the delta wing configuration (figs. 30 and 33) the total damping ratio 

will be about 0.17 higher at M = 0.8 and 0.07 higher at M = 1.2. 


14 


MCA RM H56A06 


effect (fig. 54). The reason for these trends is, of course, that the 
negative stabilizer increases the pitching velocity (fig. 55) while the 
positive input appreciably reduces the pitching rate. Reference 8 pre- 
sents flight data showing similar effects of stabilizer inputs. When the 
basic pitch damping is increased about sixfold, much less critical motions 
are experienced for the unfavorable stabilizer input (symbols, fig. 54). 

Data not presented indicated that somewhat more positive stabilizer 
inputs than shown in figure 54 were undesirable. This was attributed to 
the fact that large negative pitching rates were produced. 

It should be noted that, whereas negative stabilizer input aggra- 
vated the motions in this instance, conditions in which the basic pitching 
velocity were negative or instances in which "complementary" sideslip 
were present might be made more critical with positive stabilizer inputs 
(table I). 

It should also be pointed out that large changes in pitching moment 
due to sideslip could also have a decided effect on the characteristics. 

Effect of initial angle of attack .- The effects of initial angle of 
attack for the swept-wing configuration are presented in figure 56 for 
a principal axis inclination of 1° at two levels of directional stability. 
Data are presented at a = 5°> 1°> and -5° or at angles of attack of the 
principal axis of +4° and 0°. The results shown in figure 56 indicate 
that the peak excursions occur at about the same roll rate for negative 
and positive conditions of angle of attack with the amplitudes appearing 
as nearly mirror images of each other. This is even more evident in the 
representative time history plots of figures 57(a) and 57 (b), which illus- 
trate the expected diametric opposite nature of the motions. For this 
flight condition the major disturbing function is produced by the angle 
of attack of the principal axis through the Op term and the I^z terms 

in the equations of table I ^C np = Cng^ = oj . Thus, when the principal 

axis is lined up with the relative wind, essentially no motion is produced. 
If, however, a small yawing moment is introduced by assuming C np = -0.05> 

then relatively large peak^ amplitudes are obtained, particularly in 
regard to sideslip motions (symbols, fig. 56). This is not too surprising 
if it is considered that near p = 1.6 radians/sec an effective insta- 
bility is present for the Cn^ = 0.057 condition. Thus a relatively 

small out-of-trim movement could be expected to produce a sizable effect. 


^The yawing moment introduced would have produced about A3 = 2° 
for a roll rate of 1.6 radians if a three-degree-of -freedom uncoupled 
motion is assumed. 
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Similar trends on a somewhat reduced scale are indicated for 

c np = 0.11^ (fig. 36). 

A summary plot is presented in figure 38 for a principal axis incli- 
nation of 5°. The overall effects are, as might he expected, with the 
greatest amplitudes recorded for conditions having the largest initial 
angle of the principal axis. In figure 39 "the results of figures 3 6 and 
38 are compared hy plotting the maximum values shown for any conditions 
as functions of the initial angle of attack of the principal axis . The 
importance of the latter parameter in correlating the results is at once 
evident upon inspection of figure 39* The amplitudes were in each case 
a minimum when Op 58 0°. 

It should be mentioned that if rolling performance at initial 
g levels of 1 or greater is considered to be of primary importance, then 
a moderately large nose-down inclination of the principal axis (posi- 
tive e) would be beneficial because it would reduce otp. 


It is believed that if an appreciable value of 



or 



is 


present, the angle of attack (fig. 39) at which minimum motions occur 
would be shifted, depending on the sign and magnitude of the additional 
disturbing functions. 


The sensitivity to change in the initial angle of attack is illus- 
trated for the delta-wing configuration at a Mach number of 1.2 (fig. 4o). 
Data are presented for initial angles of attack of 2.0° and 3*9° 

(op = 0° and I. 9 0 ) • On first inspection it might seem surprising that 

such a large motion would be obtained for the case of op = 0° at high 
roll rate. However, a small value of Cnc was used (table II ) which 

provided the necessary disturbing function. When C nc was assumed 

5 h 

zero (symbols, fig. 40) , essentially no a or (3 motions were obtained. 

Effect of altitude .- The influence of a large altitude change on 
the motions developed in rolls is shown in figure ^1. Data are presented 
for the swept-wing airplane ^ Cn^ = 0 . 057 ) at initial angles of attack of 

2° and 5° and altitudes of 10,000 feet and 32,000 feet. It should be 
noted that the initial angles of attack of 2° and 5° correspond to level 
flight conditions at 10,000 feet and 32,000 feet, respectively. It is 
evident from the results that the major effect of decreasing altitude 
(at constant initial a) is to delay the buildup of large motions to a 
higher average roll velocity. From the loads standpoint, of course, 
much more critical conditions could be obtained at the lower altitude 
because of the two- and one-half -fold increase in dynamic pressure . The 
lower resonant frequency calculated from equation (l) is noted by ticks 
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for each altitude and it is apparent that the increase in average roll 
velocity at which largest coupling effects occur is predicted rather well 
from the simple theory. It might be mentioned that the decrease in maximum 
amplitude at the lower altitude is possibly attributable to the 60 percent 
increase in damping ratio in pitch assumed invariant with altitude). 

Effect of yawing moment due to roll .- For most of the calculations 
a value of Cnp = 0 was assumed for the swept-wing configuration. Fig- 
ure 42 shows the effects of large arbitrary increases in C n . for 

C no = 0.057 and 0.1l4. For the lower value of C n . an increase in C n 
P xx p 

to 0.10 is seen to reduce appreciably the size of the motions obtained. 

A further increase to 0.20 produces rather violent motions at a relatively 
low roll rate. For a value of C n p = 0.1l4 an increase in C n ^ to 

either value has a deleterious effect on the maximum amplitudes. 


It was mentioned in the discussion of initial angle-of-attack effects 

that changes in C n or C n ~ which might alleviate the motions for 

P a^ 

certain could aggravate the coupling problem for other initial angles 


O 


of attack. Figure 43 shows the effects of introducing a value of 
C n =0.10 (c n g = 0 . 057 ) in rolling maneuvers performed at cx,q = 5° 
^ o ^ o 

and -3 • For <Xq = 5 the favorable effect of increasing C n p is 


clearly shown (figs. 43(a) and (b)) and probably results from the large 
change in the yawing velocity buildup. When the roll is made at <Xq = -3°, 

however, the positive Cnp is seen to aggravate the motion appreciably. 


Effect of C ^ • - A number of Cj variations were investigated for 

H H 


the basic swept-wing configurations = 0.057)* Figure 44(a) illus- 
trates the effect of doubling the ordinates of the variation with 

a (table II ). No appreciable change in the amplitude of the motions 
is indicated. 


Unswept-wing configurations are generally characterized by much 
smaller variation in with a than are swept-wing configurations. 

It was decided, therefore, to study the effects on the roll coupling 
problem for three constant levels of C^. Figure 44(b) presents the 

results for C 7 values of 0, -O.065, and -0.126 per radian. The data 

indicate that as the level of C 7 is increased the peak motions occur 

at a slightly higher roll rate, but the peak amplitudes are not much 
reduced. Another aspect of the problem results from the combination of 
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high dihedral effect and adverse sideslip which reduces the average roll 
velocity attainable for a given aileron deflection (fig. 44). Thus, for 
the highest C 7 level if maximum aileron deflection is limited to 30° 


(although -l4° of sideslip would be obtained), the angle -of -attack excur- 
sions would be greatly reduced at the maximum obtainable roll rate of 
-I .60 radians per second (fig. 44(b)). 


Effect of mass distribution .- The importance of the mass distribu- 
tion parameter ~ treated in this section. In figure 45 data 

1 Z 

are presented illustrating the effects obtained in changing from a 

configuration with its mass concentrated along the fuselage ~ = 0.9 

i r - i x Iz 

to a condition of more equitable distribution —4= = 0.2. In this 

I Z 

study the inertia in yaw Ig was held constant (table III). 


For a C nD level of 0.057 it would appear that reducing the mass 
P I Y - Iv 

concentration along the fuselage, such that -~= — is reduced from 

iz 

0.90 to 0 . 50 , results in a slight decrease in the peak amplitude of the 
motions and an increase in the roll velocity at which these peak values 
occur. A further reduction in the inertia ratio to 0.20 appears to 
improve the situation appreciably. To compare these trends with simple 
theory, a stability plot has been prepared (fig. 46). The lower resonant 
frequency is indicated by the ticks on figure 45. As has been demon- 
strated previously in this paper, the simple theory adequately predicts 


the trends. As the mass distribution factor 


I Z 


is reduced, the 


resonant frequency is of course raised and the intersection of the sta- 
bility line is changed from one indicating an appreciable (3 divergence 
to an intersection tangent to the a divergence side of the boundary. 

A series of representative time histories showing the basic character 
of the motions for near peak coupling conditions is presented as fig- 
ure 47. Note that the large value of associated with the lowest 


inertia ratio reduces the rolling acceleration to such a degree that a 
much higher aileron deflection is required to obtain a given average 
roll rate in the 3^0° maneuvers. Otherwise, the general character of 
the motions developed is quite similar. 


For the Cnp level of 0.1l4 per radian, the summary plot (fig. 45) 

indicates very little difference in the inertia range from 0.5 to 0.9* 

I Y - Iv 

Moreover, for the condition of most equal mass distribution, — =■= — = 0.2, 
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the results are seen to he slightly worse than for the other conditions. 
The reason for this condition can he explained with the aid of the sta- 
bility chart of figure 46. For this stability level, as the mass is pro- 
gressively more evenly distributed, the characteristic stability line 
intersects further in the a divergence side of the boundary. 


General Design Considerations 

In the foregoing analysis various individual effects have been 
examined relative to the roll coupling problem at several flight condi- 
tions. A contemporary swept-wing airplane designed for Mach numbers near 
1.7 is now examined briefly to assess the roll coupling problem and also 
to determine how the results of the previous sections apply to the airplane. 

It should be noted in the airplane design under consideration (tail C 
of ref. 8) the vertical tail size was primarily dictated by static stability 
considerations at design supersonic Mach number . The longitudinal stability 
corresponded to a static margin of about 10 percent at subcritical speeds 
with the usual large increase in stability occurring at transonic speeds. 

First, a roll stability diagram is constructed (fig. 48(a)), to illus- 
trate the effect of Mach number on the stability proportioning. The deriv- 
atives used were taken from unpublished data (table TV). The stability 
diagram indicates that at M = 0.73 the proportioning of longitudinal to 
directional stability is near optimum, whereas at M = 0.93 and M = 1.26 
the normal rearward aerodynamic -center shift produces a less desirable 
intersection. The flight envelope of the airplane is presented as fig- 
ure 48(b). The solid lines represent lines of constant lower resonant 
frequency (calculated from data of table IV) and the dashed line indicates 
the maximum capabilities of the ailerons in relation to the average roll 
velocity obtainable in 360° rolls from initial lg flight. 

For this airplane, in almost the entire supersonic flight range, the 
average roll rates available are considerably lower than the roll velocity 
at which the most violent maneuvers might be expected. At moderate sub- 
sonic speeds the ailerons are capable of producing roll rates in which 
serious roll coupling could be obtained; however, the desirable propor- 
tioning of longitudinal and directional stabilities would be expected to 
have an alleviating effect on the amplitudes of the motions. 

The use of figure 48 for a quick assessment of the roll coupling 
problem is obvious. It should be realized, however, that conditions 
falling slightly to the right of the available aileron power line can also 
produce large motions as indicated in previous sections of this paper. 

This condition is particularly important when the increase in dynamic 
pressure is considered. Preparation of a figure similar to figure 48 at 
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a number of g levels supplemented by a knowledge of principal axis incli- 
nation would appear to be a logical first step in arriving at critical 
conditions for detailed analog calculations. 

As indicated in reference 8, a comprehensive roll investigation of 
the airplane (full aileron 360° rolls) has been made particularly at the 
flight conditions shown by the symbols in figure 48. At supersonic speeds, 
no serious coupling problem has been encountered in flight. At M = 0.73> 
although no violent maneuvers were experienced, the motions obtained were 
fairly large and were sensitive to inadvertent longitudinal control motions. 
At M = 0.93 fairly large motions were also obtained and calculations indi- 
cated that, had slightly higher aileron deflections been available, a rather 
serious maneuver might have been encountered. 

It would appear that the airplane of reference 8 is acceptable from 
the overall roll coupling standpoint. Based on the results of figures 26 
and 34, however, a substantial increase in pitch damping would materially 
improve the situation at subsonic speeds both for rolls made with and 
without inadvertent stabilizer inputs . At supersonic speeds the aileron 
power was considered more than adequate by the pilot, although the average 
roll rates were slightly lower than at subsonic speeds. 

To summarize, then, it would appear that present design practice 
(adequate static directional stability at design high-speed Mach number 
and moderate longitudinal stability at subcritical speeds) would probably 
furnish a good starting point in obtaining acceptable motions in rolling 
maneuvers. It might be desirable, however, to insure the presence of an 
adequate margin between the lower resonant frequency and the average roll 
velocity obtainable from the ailerons at supersonic speeds. In this con- 
nection it should be stated that the ailerons should be designed so that 
excessive roll rates are not attainable particularly at supersonic speeds. 
Finally, if a pitch damper producing a total damping ratio on the order of 
0.7 to 1.0 at subcritical speeds is used, the tendency for serious roll 
coupling effects in this speed range will be minimized. 


CONCLUSIONS 


A fairly comprehensive five-degree-of -freedom analog study of the 
roll coupling problem for a generalized swept-wing and for a tailless 
delta-wing airplane configuration has been completed and the following 
conclusions are in order: 

The various aerodynamic and inertia parameters considered indi- 
cated: 

(a) The relationship of the longitudinal and directional stabilities 
is of paramount importance. An opti m u m condition exists when the resonant 
frequencies in pitch and yaw are approximately equal. 
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(t>) Increases in pitch damping had a pronounced favorable effect in 
reducing the amplitudes of the motions encountered and were, in general, 
considerably more effective than corresponding changes in yaw damping. 

(c) Changes in mass distribution in the practical range produced 
only relatively small effects. 

(d) Alterations in the yawing moment due to rolling Cnp although 

sometimes producing favorable effects would not be useful in alleviating 
the coupling problem throughout the angle -of -attack range. 

(e) Changes in the rolling moment due to sideslip parameter C^p 
could produce significant improvements under certain conditions. 

2. A study of the effects of flight condition revealed that: 

(a) The amplitude of the motions developed depends to a large extent 
on the duration of the maneuver. Limited studies indicated that 90 1 ” 1 roll 
maneuvers would be considerably less critical than 360 polls . 

(b) ftmal 1 inadvertent stabilizer inputs can greatly affect the motions 
that would otherwise be obtained. 

(c) The engine gyroscopic terms can cause rolls in one direction to 
be somewhat more critical than corresponding rolls in the opposite direc- 
tion. 


(d) The initial angle of attack of the principal axis has an impor- 
tant bearing on the amplitude of the motions, particularly in the absence 
of other disturbing functions. 

(e) For a given initial angle of attack a reduction in altitude will 
delay critical conditions to a higher roll rate, but the maximum ampli- 
tudes may be affected only slightly. 

(f) The effects of variations in Mach number cannot easily be gener- 
alized since there are likely to be changes in most of the controlling 
parameters . 

3 . Utilizing simple concepts proved useful in assessing the quali- 
tative effects of many of the aerodyn am ic and inertia parameters and 
changes in flight condition. The calculated lower resonant frequency 
generally corresponded to the average roll velocity at which the more 
serious motions could be expected. 
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4. It is reasoned that rational design procedure can avoid the 
problem of serious roll coupling at supersonic speeds and minimize the 
problem at subsonic speeds for the configurations of the type considered. 


High-Speed Flight Station, 

National Advisory Committee for Aeronautics, 
Edwards, Calif., December 20, 1955* 
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APPENDIX A 


If the equations at the top of page 9 of NACA TN 1627 (ref. 2) are 
written to include the engine inertia terms and actual value of 

~ the equations become in the notation of reference 2 
X Y 


e - p 0 ^ - (Pot + Po 2q ) f ' + 2 ^e co ePo (® - Po^ + 

(ty + Pq®) K + O)0 2 Pq 2 9 = 0 


( 1 ) 


where 


¥ + P 0 ® + (p Q 2 t - Pq®) F + 2 ^P 0 (t + P O 0 ) - 


(0 - Pq+Jk^l + ayipQ 2 * = 0 


( 2 ) 


! X - Iy 


= F 




= K 


e 
♦ 
P r 


I* - I 


X _ 


= F' 




= K n 


X Y I Z 

angle of pitch relative to flight-path direction, radians 
angle of yaw relative to flight-path direction, radians 
steady rolling velocity, radians/sec 


Then, as in NACA TN 1627, assuming 


t ' = p n t and D = — — = — — — 

dt' p Q dt 


equations (l) and (2) become 


D 2 0 - Dilr - Di|rF' - 9F' + M + 2? 0 cs 0 D0 - 2^ e CD 0 ^ + oa 0 2 0 = 0 (3) 


o K-i D0 K-i \|/" o » 

® V + D0 + \jrF - D0F + P^yn^D^ + 2^^aY© + ny 2 ^ = 0 (4) 
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The determinant of equations (3) and (4) is 


D 2 - F' + 2^q03qD + o) fl 2 + £- 


p 0 


D _ DF . + ^ _ 2£ 0 a> 0 


K-iD 

D - DF + - U- 


+ F + 2?^ + + pj 


= 0 


The determinant may be expanded to give the quartic 

aD^ + bD^ + cD 2 + &D + e = 0 

The divergence boundary is obtained by setting the coefficient e = 0 


= -FF' + cd 0 ^ 2 - F'o)^ 2 + o) Q 2 F + + 


M Itoj, 2 KiK F'Ki Kjcoe 2 _ 
P 0 P 0 P 0 2 P 0 p 0 _ 

if £0 = = 0 then equation ( 5 ) becomes 


(5) 


= -FF' + <ne 2 afy 2 - F'city 2 + u>0 2 F + + 


KjK F'K X | _ q 

p 0 2 p 0 p 0 
from which it can be shown 


o) fi 2 = F' - = F 

9 P 0 


\ r Y p 0 


p 0 


03 2 = -F-^=-F-f ^ 


\ p zPo 


2 . I X e co e 
I Z p 0 




(6) 


(7) 


2 ly^oig 

if o)q + — is plotted as ordinate and 03 ^ + -^r— — as abscissa. 


2b 


MCA RM H 56 AO 6 


it is seen that a stability plot can be obtained similar to those in 
MCA TN 1627 


-F 




+ 


I X e co e 

!zPo 


It will be of interest to determine the roll velocity at which resonance 
occurs. The square of the undamped natural frequency in yaw can be 
written 



C np \ pV 2 Sb 


(8) 


then, substituting equation (8) into equation (7), enables calculation 
of the steady roll velocity at which instability first occurs 


? I V e 

P 0 2 + (i7Ti 


Y 


P 0 + 


Cn p \ PV 2 Sb 

TT - 


= 0 


I X e cu e 

2(I X - I Y ) 


MIX - 1 Y, 


Cnp \ PV 2 Sb N 

ix - lY 


(9) 
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Inasmuch as the engine term inside the radical is usually small relative 
to the directional stability term, equation ( 9 ) can be simplified to 


applicable when the directional divergence boundary is intersected first 
as roll velocity is progressively increased. Similarly it can be shown 


applicable when the pitch divergence boundary is first crossed. 

These equations indicate that left rolls will have a somewhat lower 
resonant frequency than right rolls for right-hand engine rotation. 

For the practical condition of finite damping in pitch and yaw, it 
appears the boundaries can be calculated as if engine terms were neglected 
(ref. 2) and then used with the modified parameters to account for the 
engine inertia effects. 



(10) 



(11) 
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APPENDIX B 


The effects of roll rate and the basic proportioning of longitudinal 
to directional stability on divergence rate can be ex am ined by treating 
a simple condition of zero damping and by neglecting the engine terms of 
equations 1 and 2 of appendix A. If these assumptions are made, the 
determinant of equations 5 and ^ of appendix A can be expanded to give 


aD 1 *- + cD 2 


+ e 


= 0 


where 


a = 1 

c = 1 - FF' + <xy 2 + odq 2 
e = -FF' + (JQq 2 ^ 2 - o^ 2 F' + <^0 2 F 
This equation can be solved as a quadratic in D 2 as follows: 



The value of the positive real root D of the preceding equation 
will determine the nondimens ional time to double amplitude in accordance 
with the formula 

t- = 0-^93 

2 Positive real root 

where t 1 ^ is the nondimens ional time to double amplitude. 
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Consider the specified points in the unstable portion of the following 
stability chart 


-F 




If -F = 0.71 and F' =0.95 then the positive real root for the 
points considered will be as follows: 


oty 2 

CD0 2 

Positive real root 

t * 2 

0.5 

2.0 

0.228 

3.03 

.5 

4.0 

.519 

2.17 

.5 

6.0 

.357 

1.94 
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Thus it would be expected that if the vertical displacement of the 
three points considered is caused, for example, by variation (the 

rate of roll therefore being identical) the higher the location of the 
point in the unstable region the greater the rate of divergence. 

If on the other hand the stability level of both and C n p 

are altered such that a given point in the unstable region is attained 
at a different roll rate, it would appear that the rate of divergence 
would be proportioned to roll rate inasmuch as 


t ' 2 = p 0 t 2 


t 2 


V 

Po 


where tg is the time to double amplitude measured in seconds. 
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TABLE I 


O 


EQUATIONS OF MOTION (REFERENCED TO BODY AXIS) 


i - * if * + if M * ft c * 8 „> * ft * ffe c v * 4fe <v * ft <v * ftfe <v 


• / Iz ” % 

<1 = pi^ ,pr + 


L X x X ^£4. 

YjCZ _2 X XZ _2 Ix e me 


i x -r 


I X ~‘P" 


i Y - i Y 


r = 


/ Ix “ ly 

(■v 1 


pq + 


X XZ • T XZ 


qr + 


I X e 0J e 


q + ^ CW 6 r + 


Ab 2 


V *V* T m. z C % r + 2VI^ V * 15 ^pP * 2VI^ ^ n pP + ^ng^at 


Ab^ 


C n J? + 4r C no |3 + C n J + C nR 8, 


Ab 


3-|m 5 -r + ap+^ C Y P + Cy P + Cy r + 


Ab 


Ab 


W "YpP - ^2 "Y P P - ^2 "Y r < - ^2 "Yp 


CyiP 


d = q + f n, - pp - A. Cla a + ^ ^ C^i* 


Zj = m^r - n^q 
mj = n^p - Z^r 
n 5 = Z^q - m 5 p 

where 

A = ^ 

A 2 


Z^ = -sin e e 

= sin cp e cos 0 e 
n^ = cos cp e cos 0 e 


cp e , 0 e angles between the body axis and earth gravity axis 
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TABLE II 


DERIVATIVES AND CONSTANTS REQUIRED FOR BASIC CALCULATIONS 



Swept -wing airplane 

Delta-wing airplane 

M = 0.7 

M = 0.8 

M = 1.2 

Basic flight condition 

hp, ft 

32,000 

40,000 

40,000 

| PV 2 , lb/ sq ft . . 

197 

175 

395 

V, ft/sec 

690 

777 

1,165 

a Q , deg 

5-0 

8.0 

3.9 

*0’ z 

1.0 

1.33 

1.0 

T 0 > deg 

0 

0 

0 

deg 

0 

0.8 

1.0 

Physical characteristics 

C , f t 

11.3 

23.13 

b, ft 

36.6 

37.8 

S, sq ft 

377 

661.5 

m, slugs 

745 


745 

Iy, slug-ft^ . . . 

10,976 

13 , 

200 

Iy, slug-ft 2 . . . 

57,100 

106,000 

I z , slug-ft 2 . . . 

64,975 

ll4,6oo 

IXZ> slug-ft 2 . . . 

942 

3,540 

e, deg 

1.0 


2.0 

I Xe t %' slug-ft 2 . . 

17,554 

17,554 
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TABIE II.- Concluded 

DERIVATIVES AND CONSTANTS REQUIRED FOR BASIC CALCULATIONS 



Swept -wing airplane 

Delta-wing airplane 

M = 0.7 

M = 0.8 

M = 1.2 

Aerodynamic derivatives for basic flight conditions 

C Li , per radian 

0 

0 

0 

Cj^j per radian 

3.88 

2.78 

1.83 

Cmi t > Per radian 

-1.0 

* 

* 

Cmq, per radian 

-3.5 

-1.5 

-0.7 

Cm^ per radian 

-0.36 

- 0.18 

-0.54 

Cj^, per radian 

-1.5 

-0.5 

-0.3 

Cmp> Per radian 

0 

0 

0 

C Q p, per radian 

0 

0 

0 

c n r > P er radian 

- 0.095 

-0.140 

-0.150 

C n p# per radian 

0.057 

0.070 

0.051 

C n i, Per radian 

p 

0 

0 

0 

C nR > Per radian 

°at 

0 

-0.014 

-0.011 

C n& ^, Per radian 

* 

* 

* 

Cy p ^ per radian 

0 

0 

0 

Cy r > per radian 

0 

0 

0 

Cyp> per radian 

- 0.28 

- 0.57 

-0.10 

Cy*, per radian 

0 

0 

0 

C per radian 

- 0.255 

- 0.195 

- 0.232 

C i T ) per radian 

0.042 

0.070 

0.080 

C , per radian 

0.054 

0.0715 

0.050 

C 7( . f per radian 

®r 

* 

* 

* 

Cj^, per radian 

0 

0 

0 

Cjp> per radian 

(See curve) 

(See curve) 

(See curve) 


*Not applicable 




a, radians 

(b) Delta-wing airplane. 
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TABLE III 


MOMENTS OF INERTIA USED IN MASS DISTRIBUTION STUDY 


I Y - I X 

x z 

Iy, slug-ft 2 

Iy, SlUg-ft 2 

I^, slug-ft 2 

0.20 

27,000 

1 + 0,000 

64,975 

• 50 

18,200 

51,000 

64,975 

•71 

10,976 

57,100 

64,975 

• 90 

4,1+30 

63,000 

64,975 


TABLE IV 


VALUES OF C n AND C m USED IN COMPUTATIONS 
0 a 


M 

S 


0.6 

0.100 

- 0.39 

• 7 

.100 

-.39 

.8 

.100 

-.43 

• 9 

.110 

-.80 

1.0 

.ii4 

- 1.14 

1.2 

.103 

-1.14 

1.4 

.063 

- 1.00 


Note: The physical and mass characteristics used to calculate 

figure L8 were the characteristics of the swept-wing 
airplane of table II. 
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r, radkms/sec 


ft deg 


a, deg 



(a) Input A; &a t = -lV 3 









Figure !•- Representative time histories 
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E-1841 

Figure 2.- Photograph showing control stick and oscilloscope presentation. 
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Figure 3.- Representative chart defining stability boundaries. 
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A/3, deg 


A a, deg 



p, radians/sec 


Figure 4.- Effect of 



variation. Swept-wing configuration; 
= 0.057} input B. 
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Figure 5.- Effect of variation. Swept-wing configuration; 

C n p = 0.114; input B. 
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Figure 6.- Effect of variation. Swept-wing configuration; 

Cnp = 0.228; input B. 



Figure 7.- Stability boundaries for swept-wing configuration. Left rolls. 
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Figure 8.- Representative time histories for swept-wing configuration. 

Input B; C n p = 0 . 057 * 
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Figure 9* _ Representative time histories for swept-ving configuration. 
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Figure 10.- Representative time histories for swept 
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-wing configuration. 
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Figure 11.- Relative importance 
Swept -wing configuration. 


of various terms in equations of motions. 
Input Bj = -0.72; C n p = 0.057* 
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Figure 12.- Relative importance of various terms in equations of motions. 
Swept-wing configurations. Input B; = -0.l8j C n ^ = 0.228. 
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Cj^ variation. Delta-wing configuration; 
M = 0 . 8 ; input A. 


Figure 13.- Effect of 
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Figure lA.- Effect of and C n p combinations. Delta -wing 

configuration M = 1.2; input A. 
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Figure 15.- Stability boundaries for delta -wing configurations. M = 0.80. 

Left rolls. 
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Figure l6. - Stability boundary for delta -wing configuration. M = 1 . 2 . 
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(a) Input B; C n p = 0.057* 
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Figure 17.- Effect of roll direction. Swept-wing configuration 
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Figure 19.- Effect of roll duration. Swept-wing configuration. 
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Figure 20.- Representative time histories showing effects of roll duration. Input B; C n p = 0.057* 
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Figure 21.- Representative time histories showing roll duration effects 
Delta -wing configuration. Input A; M = 1.2. 
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Figure 22.- Effect of increasing C nr . Swept -wing configuration. 

Input B; C n p = 0.057* 
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Figure 23.- Stability chart illustrating effects of increasing C nr . 
Swept -wing configuration. Left rolls. 
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(b) c nr = - 0 . 95 . S = -9-5- 


Figure 2k, - Time histories showing effect of increasing C nr . Swept 
wing configuration. Input Bj = -0.72 j C n p = 0.057* 
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(a) C nr = -0.095. 





Figure 25.- Time histories showing effects of increasing C nr . Swept- 
wing configuration. Input Bj = -0.18; C n p = 0.228. 
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Figure 26.- Effects of increasing C m . Swept -wing configuration. 
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(a) C np = 0.057. (b) C np = 0.114. 


Figure 27 . - Stability charts illustrating effects of increasing C m . Swept-wing 

configuration. Left rolls. 
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Figure 28.- Time histories showing effects of increasing C m ^. Swept- 
wing configuration. Input B; = -0.72j C n p = 0.057* 
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Figure 29 .- Time histories showing effects of increasing C m q. Swept- 
wing configuration. Input Bj Cn^ = -0.l8j C n p = 0.228. 
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Figure 30-“ Effect of increasing C^. Delta -wing configuration. 

M = 0.80; input A. 



Figure 31.- Stability chart illustrating effects of increasing C^. 
Delta -wing configuration. M = 0 . 80 . Left rolls. 
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Figure 32.- Effects of increasing C m ^. Delta -wing configuration. 

M = 1.20; input A. 



Figure 55. - Stability chart illustrating effects of increasing C m 
Delta -wing configuration. M = 1.20. Left rolls. 
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Figure 34.- Effects of simultaneous aileron and stabilizer inputs. 
Swept-wing configuration. C n p = 0.1l4j input A. 
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Figure 35.- Time histories illustrating effects of simultaneous aileron and stabilizer inputs. 

Swept-wing configuration. Input A; C n p = 0. 114. 
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(b) Input A; C np = O.llU. 


Figure 56 .- Effects of varying initial angle of attack. Swept-wing configuration, e = 1° 
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Figure 57* - Time histories showing effects of different initial angles 
of attack. Swept-wing configuration. e = 1°. 
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Figure 58 .- Effects of varying initial angle of attack. Swept-ving configuration, e = 5 ° 
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Figure 39 * - Variation of peak motions with initial angle of attack of 
principal axis. Swept-wing configuration. 
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Figure AO.- Effects of initial angle of attack. 

tion. M = 1 . 2 . 
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Figure 4l.- Effect of altitude. 
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Figure h-2.- Effects of increasing 




-4 -.8 -1.2 -1.6 -2.0 -2.4 -2.8 -3.2 

p, radians/sec 


(b) Input A; C n p = 0.114. 


'V 


Swept-wing configuration 




NACA RM H 56 AO 6 


Figure kj 


Sq| , deg 



£,deg 


a ,deg 












- 20 - 




C\ 






(a) Qq = 5°; = 0. (b) Oq = 5°; * 0.10. 
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Time histories showing effects of C n increase at two values of initial angle of 

hr 

attack. Swept-wing configuration. C n ^ = O.O57. 
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Figure 4-5.- Effects of variations in mass distribution 
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Figure k6.- Stability charts shoving effects of variations in mass distribution. Swept-wing 

configuration. Left rolls. 


— 3 
VO 


NA.CA RM H56A06 


80 


NA.CA RM H56A06 


S Qf ,deg 



r, radians/sec 0~ — — 4= 


, 0 U \ Vi t m 


/3,deg 


a, deg 










<*> h ^ L ‘°- 9 °- 


Figure 47- - Time histories showing effects of variations in mass distri- 
bution. Swept-wing configuration. C n p = 0.057* 
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Figure lj-8. - Preliminary roll analysis. Swept-wing configuration. Initial level flight. 
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